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We present a search for a new narrow, spin-1, high mass resonance decaying to fi'^fJ,~ + X, using 
a matrix element based hkelihood and a simultaneous measurement of the resonance mass and 
production rate. In data with 4.6 fb"'^ of integrated luminosity collected by the CDF detector in 
pp collisions at ^/s = 1960 GcV, the most likely signal cross section is consistent with zero at 16% 
confidence level. We therefore do not observe evidence for a high mass resonance, and place limits 
on models predicting spin-1 resonances, including M > 1071 GeV/c^ at 95% confidence level for a 
Z' boson with the same couplings to fermions as the Z boson. 
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dieted generieally in models with additional gauge groups 
a feature of many extensions to the standard model 
of particle physics (SM). Typical examples are Little 
Higgs models [2| and the next-to-minimal supersymmet- 
ric model Q. 

Current 95% confidence level (CL) lower limits on the 
mass of a Z' boson with the same couplings to fermions 
as the Z boson iZ'gj^j) are 1030 GeV/c^ from a search 
of CDF dimuon data with 2.3 fb~^ of integrated lumi- 
nosity 0, and 1023 GeV/c^ fH^ and 963 GeV/c^ [6] from 
respective searches in dielectron data from DO in 5.4 fb^^ 
and from CDF in 2.5 fb~^, respectively. 

This Letter reports a new search of the CDF dimuon 
data with several significant enhancements: twice the in- 
tegrated luminosity, a matrix element based likelihood 
providing an approximately 20% relative increase in cross 
section sensitivity at large Z' mass, and a new statistical 
approach designed to maximize simultaneously both dis- 
covery potential and mass exclusion limits in searches for 
new physics. This new approach directly fits the data to 
a single Z' boson hypothesis in the plane of boson mass 
and signal cross section, and can be applied to any search 
for a hypothetical new particle of unknown mass. 

We use a data sample corresponding to an integrated 
luminosity of 4.6 fb~^, collected with the CDF II detec- 
tor , a general purpose detector designed to study 
collisions at the Fermilab Tevatron. The tracking system 
consists of a cylindrical open-cell drift chamber and sili- 
con microstrip detectors in a 1.4 T magnetic field parallel 
to the beam axis. The silicon detectorsprovide tracking 
information for pseudorapidity jTyj < 2 [8| and are used to 
reconstruct collision and decay points. The drift chamber 
surrounds the silicon detectors and gives full coverage in 
the central pseudorapidity region (1 77 1 < 1). For the muon 
kinematics relevant to this search, the tracker provides an 
invariant mass resolution of 5M/M'^ « 0.17%/TeV/c2. 
Electromagnetic and hadronic calorimeters surrounding 
the tracking system measure energy from particle show- 
ers and minimum ionizing particles such as muons. Drift 
chambers and scintillators located outside the calorime- 
ters detect muons in the central pseudorapidity region 

H < 1- 

Events used for this search are selected online with 
a trigger requiring a muon candidate with px > 18 
GeV/c. The event selection is unchanged from the previ- 
ous search [3], requiring at least two oppositely charged 
muons with pT > 30 GeV/c and no identified cosmic 
rays Each muon is required to have calorimeter de- 
posits consistent with a minimum ionizing particle and 
have a track fully within the fiducial volume of the drift 
chamber. At least one muon must have an associated 
muon-chamber track. 

The dominant and irreducible background is standard 
model Z/j* production with subsequent decay to muons. 
We model this background using events generated at 
leading order by pythia [i3| with the CTEQ6l 11 1 par- 



ton distribution functions (PDF). Events are reweighted 
by the ratio of next-to-next-to-leading order to leading- 



order cross sections as a function of dimuon mass [12 1. 
We use the standard CDF simulation [3] to model the 
detector response to the particles in the event. The over- 
all normalization of this background is derived from data 
near the Z boson pole mass, 70 < Af^^ < 110 GeV/c^. 

The remaining background contributions are small rel- 
ative to standard model sources of dimuons. Decays of 
WW (1%) and ti (1%) to dimuons are described using 
PYTHIA and the CDF detector simulation, normalized to 
NLO cross sections [l3|- The background from objects 
misidentified as muons (0.5%) is estimated using the dis- 
tribution of calorimeter energy in the vicinity of each 
candidate muon in the sample [l^. Cosmic ray back- 
grounds (< 0.01%) are modeled using identified cosmic 
ray events, weighted by the probability to survive the 
cosmic ray removal algorithm 

The observed spectrum of dimuon invariant mass 
{Mfifj,) is in good agreement with the total SM plus 
cosmic-ray prediction, as shown in Fig. [T] We expect 
1851 ± 90 events with M^^ > 130 GeV/c^, and observe 
1813 events. 

We model the production and decay of a spin-1 res- 
onance Z'gj^^j using MADEVENT [3], assuming fermionic 
couplings equal to those of the Z boson. Initial-state 
QCD radiation, hadronization, and showering are mod- 
elled with PYTHIA. The acceptance for Z' boson events 
varies with mass, increasing from 20% at 200 GeV/c^ 
to nearly 40% at 1 TeV/c^, then dropping for AIz' > 1 
TeV/c^ due to lack of initial-state partons with sufficient 
momentum. 

The previous CDF search [4] used a binned likelihood 
fit with bins uniform in the inverse of the reconstructed 
dimuon mass to extract the best-fit signal cross section 
at a range of masses. That approach weighted events in 
the same reconstructed inverse mass bin equally, indepen- 
dent of expected track resolutions. We improve on this 
method by using an unbinned likelihood that includes a 
theoretical model of the full kinematics of the event and 
the event-by-event knowledge of the muon px resolution. 
Well- measured events have narrower likelihoods, reflect- 
ing the higher quality of their information and making 
a stronger impact on the fitted result. The likelihood is 
then used both to set limits on the Z' cross section as 
a function of mass and to construct a 2D interval that 
gives a well-defined discovery condition. The likelihood 
is given by 

L{Mz',sz') = ^ "'^^ ^ Y[Lix,\Mz',sz'), 

i 

where Mz' is the Z' pole mass, N is the number of events 
in the data, N(,g,z' are the number of events due to back- 
ground and Z' signal, Xi represents the observed kine- 
matics of the i-th event, and sz' — Nz'/{Ni,g + Nz')- 
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FIG. 1: Top: The distributions of M^^ for data with 4.6 fb~^ 
of integrated luminosity (triangles) and expected SM back- 
grounds (histograms) with two example signals. Bot- 
tom: The relative difference between observed and expected 
data, as a function of dimuon mass. Error bars show statisti- 
cal uncertainty. 



The dependence of the per-event likelihood is given by 



L{xi\Mz',sz') = sz'Lz'{xi\Mz') + (1 - sz')Lz/~f'{xi), 

and the likelihood Lz' is calculated by integrating the 
matrix element for Z' production convolved with PDFs 
and the detector resolution functions: 



d<^{qi ,q2)\Mz'{qi,q2, Mz' ) l^fpDFfpDF 



X X T{pi,qi,(Tp^jT{p2,q2,(JpT2) x Ppriqi +92,A'jcts), 

where pi^2 represent the four- vectors of the two measured 
muons, qi^2 represent the unknown four-vectors of the 
two true muons, $ represents phase space for the true 
muons, A4 is the matrix element, fpoF is the parton 



distribution function, T(p, q, ap^) is the transfer function 
that parametrizes the detector resolution as a function of 
the measured uncertainty (Jp^ , and PpT is the probability 
density function for px of the /x/x system, parameterized 
in the number of jets (Njets) with Et > 15 GeV and 
1 77 1 < 2.5. The distributions of o-p^ for the Z' signal and 
dominant Z/j* background are the same in the phase 
space region near the hypothesized Z' mass, thus they 
do not affect the likelihood ratio ordering. The distri- 
bution of PpT is obtained from simulated samples with 
initial- and final-state radiation. An analogous expres- 
sion is used for Lz/j*, which describes the likelihood for 
the dominant Z/j* background. 

We analyze the resulting likelihood in two ways. First, 
we aim to discover the regions in {Mz',sz') that are con- 
sistent with the data and inconsistent with the SM, mak- 
ing no assumptions about the relationship between Mz' 
and sz'- We refer to this as the 2D interval analysis. 
Second, we wish to set limits on the Z' mass in specific 
models. In that case, we perform a raster scan, in which 
we choose a set of values of Mz' , and at each point derive 
limits on sz'- Together with a prediction for sz'iMz') 
in a specific theory, we can use the raster scan to place 
lower limits on Mz' ■ 

The 2D interval is constructed via the unified ordering 
scheme [l7j in two dimensions, resonance mass and cross 
section (see Fig. [2|). At each test point in the {Mz' , sz') 
space, we calculate the ratio of the likelihood at the 
test point to the likelihood at the best-fit point where 
the likelihood is maximized. To determine which test 
points are consistent with the data at a given confidence 
level, we perform pseudo-experiments to determine how 
often we expect to observe such a likelihood ratio. The 
pseudo-experiments include all backgrounds and interfer- 
ence effects between the Z and Z', as well as variation 
of the nuisance parameters from systematic uncertain- 
ties described below. This approach is well designed for 
discovery, as it tests the background hypothesis exactly 
once. The significance of an observation corresponds to 
the first confidence-level contour that includes a signal 
rate of zero. It also provides a summary of Z' mod- 
els consistent with the data without relying on specific 
model details. In the 2D interval analysis, the best-fit 
signal cross section of cr = 26 fb occurs at a resonance 
mass of M = 199 GeV/c^, but is consistent with zero at 
16% confidence level. 

The raster scan is the traditional approach used in 
many analyses, including the previous Z' search [1]. It 
provides mass limits on theories that enforce a relation- 
ship between the signal fraction and the resonance mass, 
and is appropriate if outside information indicates a par- 
ticular mass is interesting. In the presence of a significant 
excess above the background-only hyphothesis, one must 
account for the number of possible Z' masses, each of 
which tests the background-only hypothesis: the look- 
elsewhere effect. The 2D analysis needs no such correc- 
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FIG. 2: Observed 95% CL limits in ID raster scan (top) and 
2D intervals (bottom) for data with 4.6 fb~^ of integrated lu- 
minosity. The solid circle indicates the best-fit value and the 
lines define regions in which 25, 50, 68, 95, and 99% of exper- 
iments would yield results less consistent with the standard 
model and the data. 



tion. While the resulting discovery significance of the 
corrected raster scan will be correct, one will still be left 
with an interval in the signal fraction at every mass. In 
contrast, in the presence of a signal, the 2D analysis will 
provide a range of masses that is consistent with the sig- 
nal. 

TABLE I: Mass limits on specific spin-1 Z' models fT3| in data 
with 4.6 fb~^ of integrated luminosity at 95% confidence level. 

Model Z[ ^sec Z'j^ Z'^ Z'^ Z'^ Z'sm 

Mass Limit (GeV/c^) 817 858 900 917 930 938 1071 



Dominant systematic uncertainties 4] include uncer- 
tainties on the PDFs and the dependence of the next- 
to-leading order cross section on the dimuon invariant 
mass. These weaken the final limits by 5-10% depending 
on mass. Additional uncertainties are the level of initial 
state radiation and muon acceptance at large transverse 
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FIG. 3: Observed 95% CL limits for data with 4.6 fb"^ 
of integrated luminosity expressed as limits on the up and 
down type charges and Cd The solid and dotted lines 
show possible models in U(X)b-xl and [/(l)io+a;5 groups, 
respectively. The dashed lines show the range for models in 
the \J(X)q-^xu group. 



momentum. 

The raster scan in mass allows us to set strong limits 
on specific models of Z' production; see Fig. [2] and Table 
m The production cross section times branching fraction 
to the dimuon final state is determined by the couplings 
of the fermions to the Z' . Figure [3] shows how mass limits 
depend on the charges of the up- and down- type fermions 
to the U(l) group associated with the Z' . Table H] shows 
the limits for the specific models described in Ref [l2| . 

In conclusion, we have applied the matrix-element- 
based likelihood technique to a search for new spin-1 res- 
onances decaying to muon pairs, set the strongest limits 
to date on the resonance cross section and mass, and in- 
troduced a statistical analysis approach that is useful for 
this analysis as well as for potential LHC discoveries. 
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